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Abstract 

The Campbell theorem, relating the variance of the cnrrent of a fission chamber (a “hltered 
Poisson process”) to the intensity of the detection events and to the detector pulse shape, 
becomes invalid when the neutrons generating the fission chamber current are not indepen¬ 
dent. Recently a formalism was developed [1] by which the variance of the detector current 
could be calculated for detecting neutrons in a subcritical multiplying system, where the 
detection events are obviously not independent. In the present paper, the previous formal¬ 
ism, which only accounted for prompt neutrons, is generalised to account also for delayed 
neutrons. A rigorous probabilistic analysis of the detector current was performed by us¬ 
ing the same simple, but realistic detector model as in the previous work. The results of 
the present analysis made it possible to determine the bias of the traditional Campbelling 
techniques both qualitatively and quantitatively. The results show that the variance still 
remains proportional to the detection intensity, and is thus suitable for the monitoring of 
the mean flux, but the calibration factor between the variance and the detection intensity is 
an involved function of the detector pulse shape and the subcritical reactivity of the system, 
which diverges for critical systems. 


1. Introduction 

The main application area of hssion chambers is the measurement of the neutron flux in 
operating (critical) reactors [2]. Fission chambers offer several advantages: they are robust; 
can be operated in both pulse and current mode and they endure high temperatures. 

One special advantage of fission chambers is their capability of suppressing unwanted 
minority components in the detector current, such as gamma events, with a proper signal 
processing technique. This is based on the so-called Campbell theorem, which establishes 
relationships between the various order moments of the signal [31 i] , hence make it possible 
to determine the mean detection rate from the second moment of the detector current [5]. 

However, using fission chambers in Campbelling mode in measurements in a reactor has 
been controversial right from the beginning. Namely, the Campbell theorems are only valid 
in the case when the detection events are independent and the detection intervals obey an 
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exponential distribution. In the case of a multiplying system, either critical or subcriti- 
cal, the detections will not be independent, due to the branching character of the neutron 
multiplication process. In fact, it is the deviation of the detection times from the exponen¬ 
tial distribution, or the deviation of the number of detections from a Poisson distribution, 
which is utilised in the reactivity measurement techniques based on pulse counting, i.e. the 
Feynman- and Rossi-alpha techniques [6]. 

Even if it is surmised that application of the Campbell theorem might be allowable even 
for detections of time-correlated neutrons, by correcting for the (presumably small) quanti¬ 
tative error by a calibration procedure, it is of fundamental importance to understand the 
qualitative and quantitative effect of the existence of the correlations between the detection 
events on the variance of the signal. This would quantify the bias of the application of the 
traditional Campbell relationships to extract the mean detection intensity from the variance 
of the detector currents, and eventually even make it possible to extract information about 
the system, such as the subcritical reactivity, the same way as it is done by the pulse counting 
techniques. 

A hrst step to achieve these goals was made recently by the present authors by setting up 
a formalism which unites the stochastic description of the branching process with that of the 
statistical theory of the detector signal [1]. In that work, delayed neutrons were not accounted 
for. In the present work we extend the formalism to include also delayed neutrons. Naturally, 
the formalism becomes more involved, and both the derivations and the results become less 
transparent. Hence, out of the two goals of the previous work (calculating the bias of the 
traditional Campbell theorem, and calculating the auto-covariance of the detector signal in 
order to determine the subcritical reactivity similarly as in a Rossi-alpha measurement), only 
the hrst will be aimed at in this paper; the second would lead to prohibitively complicated 
expressions. 

One beneht of the more complicated calculations is that quantitative results can be 
presented about the bias of the traditional Campelling method in terms of the subcritical 
reactivity of the system, with direct relevance to operating reactors. The results show that 
the bias of the traditional Campbelling technique is a monotonic function of the subcritical 
reactivity. It increases when approaching criticality and diverges in a critical reactor. On 
the other hand, in deep subcritical systems, the bias of the traditional Campbelling method 
vanishes, and the traditional formula becomes exact in the limit of a purely absorbing system 
where no branching takes place and hence the individual detection events are independent. 

2. Basic considerations 

The formalism used in this work was elaborated in two previous publications. The 
statistical theory of the hssion detector signal due to independent detection events, based 
on the backward master equation approach, was introduced in [7]. Then the theory was 
extended to the case of detecting neutrons in a subcritical medium driven with an extraneous 
neutron source with Poisson statistics pQ . All hssion neutrons were considered as prompt in 
this latter work. Although we stall try to make this paper self-contained as much as possible, 
reference will be made to the above publications for the details whenever it is practical. 

As in the previous works, a basic quantity used will be the probability 

V{y <'ns{t)<y + dy \r]s{0) = 0} = h{y, t) dy + o{dy) (1) 


2 


that at a time instant t > 0 the valne of a single detector pnlse ?7s(t), initiated by the 
detection of a nentron at time t = 0, is in the interval [y, y + dy). We assnme that the cnrrent 
pnlse generated by a nentron arriving to the detector can be considered as a response fnnction 
of the detector. In many cases, this response fnnction cannot be given by a deterministic 
fnnction /(t); rather, it shonld be described by a fnnction t) which depends on a properly 
selected random variable Hence the fnnction represents the cnrrent signal which 

exists in the detector at time t > 0 after the arrival of one nentron at t = 0. We assnme that 
this signal depends also on a random variable which is dehned by its distribntion fnnction 
w{x) 

'P < x} = / w{x') dx', (2) 

J — OO 

thns one can write that 

^+oo 

h{y,t)= d[y - ip{x,t)] w{x) dx, y>0. (3) 


The continnonsly arriving nentrons generate the detector cnrrent as the aggregate of snch 
cnrrent signals, each related to different realizations of 
In fnrther calcnlations the characteristic fnnction 


h{u;, t) = 


b+oo 


h{y,t) dy 


(4) 


will be often nsed, which is obtained from the above as 

/ +00 / ^+oo \ r+oo 

/ / 5[y — ip{x,t)] w{x) dx\ dy= w{x) dx. (5) 

■OO \j —OO / J —OO 

The moments of r/s(f) are given by the formnlae: 



k = l,2,... . 

For any given signal shape (p{x,t) and amplitnde distribntion w{x), these moments can be 
calculated. 

The main objective of the present work is to determine the stochastic properties of the 
detector signal for detections in a multiplying system driven by a Poisson-like neutron source 
with constant intensity, where the individual detection events are not independent. In the 
previous work pQ it was already shown that the Campbell theorem becomes invalid when the 
detected neutrons lose their independence. In the present work, the theory will be extended 
for the case when the effect of the delayed neutrons in the neutron multiplication process is 
taken into account. For the sake of simplicity, only one type of precursors will be considered 
in the calculations. 

Denote by 

Viy^ riit) <y + dy,t\0} = P {y, 1 10) dy (7) 
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the probability that in a subcritical system which is driven by a Poisson-like neutron source 
with constant intensity sq, at the time moment t > 0 , the detector current 77 (t) is found in 
the interval {y, y + dy\^ provided that at the time instant t = 0 the detector current and 
the numbers of neutrons and precursors were zero. For further calculations one needs the 
characteristic function 

/ -l-oo 

e-yp{y,m dy, ( 8 ) 

■00 

which will play the role of the generating function. 

Adding up for the mutually exclusive events that there will be or will not be a hrst 
collision between [ 0 ,f] and applying the convolution theorem for the latter case, one obtains 
the following backward equation for P {y,t\0): 


P{y,t\0) = e-^°^6{y) + 


So y e jj p{yi,t'\n{0) = l)P{y2,t'\0) dyidy 2 , (9) 

yi+y 2 =y 

where 

p{y,t\n{0) = 1 ) dy 

is the probability that in a subcritical system without a source, at the time instant t > 0 , 
the detector current 77(f) lies within the interval ( 7 /, y + dy], provided that at time f = 0 the 
detector current and the number of precursors were zero, while the number of the neutrons 
was equal to 1. One can callp( 7 /, 1 177 ( 0 ) = 1) the single-particle induced distribution, whereas 
P{y,t |0) is the source-induced distribution. 

Introducing the characteristic function 


b+oo 


g{uj,t\n{0) = 1 ) = / p{y,t\n{0) = 1 ), 

J —00 

and the convolution theorem, from equation (|^ one obtains the following equation: 

Gioj, t) = + So f f' | 77 ( 0 ) = 1 ) G{uj, t') dt', 

Jo 

It is easy to show that the solution of the integral equation (0 is given by 


( 10 ) 


( 11 ) 


G{u:, t) = exp < sq 


|77(0) = 1) - 1 


Uo 


dt'\. 


( 12 ) 


It can be seen that the density function p{ii,t\nJS) = 1 ), or its characteristic function 
g{uj,t' | 77 ( 0 ) = 1 ), plays a fundamental role in the description of random behaviour of the 
detector current. 

It remains to derive an equation for p{y,t\n{0) = 1), and it is at this point that our 
description will deviate from the previous work. First of all, since both the detector signal, 
as well as the number distribution of the neutrons, will be the result of a branching process, 
one can only derive a backward master equation if one also keeps track of the time evolution 
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of the number of neutrons. In addition, we need to take into account that the branching 
will generate not only neutrons but also delayed neutron precursors, as well as that a singe 
delayed neutron precursor too, can initiate a branching process and a corresponding detector 
signal evolution. 

Hence, in deriving the corresponding master equations, we need to consider the extended 
densities 


Px{y,n,c,t\n{0) = 1 ) 

( 13 ) 

Px{y,n,c,t\c{0) = 1 ). 

( 14 ) 


Here, Pxhj) c, t |n(0) = 1) dy is the probability that in a (source-free) subcritical system at 
the time instant t > 0 , the detector current r]{t) will he within the interval {y, y + dy], while 
the number of neutrons n{t) and that of precursors c{t) are equal to n and c, respectively, 
provided that at the time instant f = 0 the detector current and the number of precursors 
was zero, and the number of the neutrons was equal to 1. Likewise, PxiUi c, t |c(0) = 1) dy 
is the same probability, except that at time f = 0 the number of neutrons was equal to 
zero, and the number of precursors was 1. Once these quantities are determined, the density 
p{y,t |n( 0 ) = 1 ) appearing in ([^ is obtained as 

OO OO 

p {y, t |n(0) = 1) = Px{y, n, c, t |n(0) = 1) (15) 

n=0 c=0 

and similarly for p (?/, t |c(0) = 1 ). 

The derivation of the backward equation for Px{y, n, c, t |n(0) = 1) goes as follows. If at 
f = 0 one single neutron exists in the system, then in the time interval ( 0 , t] four mutually 
exclusive events can take place: 

• the neutron will not have any reaction; 

• the neutron gets detected in the detector with an intensity and creates a current 
pulse; 

• the neutron is captured in the subcritical medium with intensity Ac, 

• the neutron creates a hssion in the subcritical medium with intensity A/. 

It is clear that the total intensity of a reaction in the system is A^ = A^ -|- Ac -l- Aj. Similarly 
to previous work, the fact that the detection itself is a hssion event, which also will produce 
further neutrons, will be neglected. 

In a hssion reaction, k > 0 neutrons and £ > 0 precursors of the same type are produced 
with probability f{k,£). It is assumed that the number of neutrons and that of precursors 
are independent, i.e. 

( 16 ) 

For later use, introduce the generating functions 

OO OO 

= = (17) 

k=o e=o 


5 


By applying the backward approach, one can write 


Pa:{y,n,c,t\n{0) = 1) = e ^(i/) ^c,o + Ad ^n,o 4,o / e h{y,t') dt'+ 


Ac5n,o5c,o / 




X 


where 




"’EE /<»'(4)/<“>((!) 

II ^ ^ 

k £ yi+y2=y ”i+”2=ri ci+c2=i. 

ni, Cl, t' |n(0) = 1) V^(|/ 2 , n 2 , C 2 , t' |c(0) = l)dyi dy 2 dt', 

Uk{yi, rii, Cl, t' |n(0) = !) = [!- A{k)] 6{yi) 5^,0 5ci,o+ 

■ n 1^(0) = 1) 


(18) 


(19) 


yn+-+yifc=?;i cii+-cife=ci i=l 


and 


A(£) 


Ve{y2, n2, C2, t' |c( 0 ) = !) = [!- A{£)] 5(1/2) 5 n 2 ,o <^c2,o+ 

£ 

j ■ j Y n n2j, C2j, t' |c( 0 ) = 1) dy 2 j. ( 20 ) 


y21 + -+y2l=y2 ” 21 +-" 2^="2 ^21+-C2l=C2 j=l 


In a similar manner, taking into account the two mutually exclusive events that the delayed 
neutron precursor will not decay or will decay with intensity A, the following equation can 
be derived for the case when the branching process is started by one precursor: 


P: 


,{y,n,c,t\c{0) = 1) = e S{y) 6n,o5c,i + X / e px{y,n, c,t'\n{0) = 1) dt', (21) 


which connects the density function Px{y, n, c, t |c(0) = 1) with Px{y, n, c, t |n(0) = 1). 
Dehning the generating functions 


gx{uj,zi,Z2,t\n{0) = 1 ) = 


/ + CXD 

^zujy Y p^{y^ n, c, t |n(0) = 1) z'^ z^ dy (22) 

00 ^_ r\ ^ _n 


and 


gx{u,Zi,Z2,t\c{0) = 1 ) = 


n=0 c=0 


/ _1_C)0 'AO 00 

^zujy Y Y Px{y, n, c, t |c(0) = 1) x” X 2 dy, (23) 

•00 ^_n-n 


71=0 C=0 


from (18) one obtains the equations for the generating functions in the following form: 

rt rt 


gx{uj, Zi, Z2,t\n{0) = 1) = e '^’’*Xi + Ad / e h{u,t') dt' + Xc / e dt'+ 

Jo Jo 

Xf [ e"^Ai-t')g(p) [g^(uj,Zi,Z 2 ,t' \n{0) = 1)] [gx{uj, Zi, Z 2 ,t' \c{0) = 1)], (24) 
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while the generating function of equation (21) is given by 


ga,{uj,zi,Z 2 .,t\c{Q) = 1) = e ^^Z 2 + X / e gx{u:, zi, Z 2 .,t'\n{Q) = 1) dt'. (25) 

Jo 

At this point it is possible to revert to the case when only the stochastic behaviour of the 
detector current is of interest only, irrespective of the number of the neutrons or precursors 
in the system. The dynamics of the branching is compressed into the non-linear functions 
[... ] and [.. 

substituting zi = Z 2 = I and by turning to the quantities 


q''^' [... ] and q^'^’ [•••]• Hence in the continuation we will simplify Eqs (24) and (25) by 


and 


g{u, t |n(0) = 1) = gx{uj, 1,1, t |n(0) = 1) 
g{u, t |c(0) = 1) = gx{uj, 1,1, t |c(0) = 1) 


(26) 


(27) 


As is seen from (12) it is the g{ui,t |n(0) = 1) of (26) which is needed for the calculation of 


from which the moments of the stationary detector current in a subcritical system 
driven by an external neutron noise can be determined. 


From equations (24) and (25), applying the notations dehned in (26) and (27), one obtains 

g{uj,t\n{0) = 1) = + Xd [ h(a;, t') dt' + Xc [ dt'+ 

Jo Jo 




(28) 


and 


g{u, t |c(0) = 1) = e -I- A / e g{u, t' |n(0) = 1) dt'. 


(29) 


For the determination of the cumulants of the detector current we will use the well-known 
relation 

T\" ra’"K(a;,f)' 


where 


f^n{t) = - 


K(a;,f) = lnG'(a;,f) = sq 


aw" 


ij=0 


g{u,t\n{0) = 1) - 1 


uo 


dt'. 


(30) 

(31) 


3. Expectation of the detector current 


By using expressions (30) and (31), one can write 

^dg{u,t' |n(0) = 1) 


1 


{r]{t)) = Ki{t) = Ii{t) = So - 


>0 L 


where 


dg{uj, t |n(0) = 1) 
duj 


du 


= ii{t |n(0) = 1) 


dt' ^ 


u )=0 


- Cd = 0 


(32) 


(33) 
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is the expectation of the detector current generated by a single starting neutron. From (28) 
one can derive the equation 


|n(0) = l) = \d / <Fi(F) dt'+ 






q[^'^ ii{t' |n(0) = 1) + ii{t' |c(0) = 1) dt\ 


id). 


where 


and 


d>i(F) = - 
% 


dh{u, t') 
du 


(34) 

(35) 


ip) 

Qi = 


dq^P\'. 

dz 


= r', 


J ct;=0 


while qi^^^ = 


pi 


J z=l 


dq^^\^ 

dz 


= ^d- 


2 = 1 


In order to obtain the solution of (34), one has to take into account the relation 

ii{t |c(0) = 1) = A [ ii{t' |n(0) = 1) dt' 

Jo 


(36) 


which follows from (29), and apply the Laplace transforms 

rt 


>0 


and 


as well as 


-st 


^—st 


^—st 


ii(t n(0) = 1) dt 

(37) 

ii{t c(0) = 1) dt, 

(38) 

4>i(f) dt. 

(39) 


It is seen that the Laplace transform of (34) satishes the equation: 


ii(s |n(0) = 1) = 


"5 + Afi + Ac + A 


^>i(s) + 


A) 


/ 


■s + A(i + Ac + A 


/ 


z/oii(s |n(0) = 1) + z/rfii(s |c(0) = 1) , (40) 


while the Laplace transform of (36) obeys the equation 

zi(s |c(0) = 1) = —^ ii(s |n(0) = 1). 

s + A 


(41) 


After elementary algebra, from Eqs (40) and (41), one obtains the Laplace transform of the 
expectation of the detector current generated by a single starting neutron in the following 
form: _ 

Ad (s + A) 4*1(5) 


ii(5 |n(0) = 1) = 


(s + A) [s + Ad + Ac + Xf{l — uq)] — X Xf Ud 


(42) 


























By using conventional notations, one has 


ii(s |n( 0 ) = 1 ) = 


Ad (s + A) *hi(s) 


(s + A) ( s + 


where 


P = 


13- P 
A 


^ A/ — (Ad + Ac + A/) 

u Xf 




(43) 


(44) 


is the reactivity, while u = Up + J^d and {3 = ^. Further, the following notations will also be 


used: 


A = 


z/A/ 


which is the prompt neutron generation time and 

13- p 


a = 


A 


> 0 


(45) 


is the prompt neutron decay constant. 

In terms of the negative roots si and S 2 of the characteristic equation 

+ s (A + a) - ^ = 0, 


expression (43) can be rewritten in the following form: 

Ad (s + A) *hi(s) 


ii(s |n( 0 ) = 1 ) = 


(s + Si) (s + S 2 ) 


where 


and 


Si = 


S2 = 


A + a + ^/ (A + a)^ + 4 ^ 


A “h Q 


(\ + af + 4^ 


(46) 


(47) 


(48) 


The expectation of the detector current, generated by a chain of neutrons generated in 
a subcritical multiplying assembly, started by one single source neutron, can be obtained by 


the inverse Laplace transformation of (|46|). It is easy to show that 

rt 


ii(t |n( 0 ) = 1 ) = Ad 


(si - A) _ y) ^-S2(t-u) 


'0 


Si - S2 


<Fi(m) du. 


(49) 


In order to obtain the expectation of the detector current in a subcritical multiplying medium 
driven by a stationary Poisson source of neutrons with intensity Sq , one has to calculate the 
integral 


= So / ii(t' |n( 0 ) = 1) dt\ 


(50) 
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the Laplace transform of which is given by the formula 


Ji(s) — So 


Ad (s + A) *hi(s) 
s (s + Si) (s + S 2 ) 


( 51 ) 


By using the Tauber theorem, one can determine the asymptotically stationary expectation 
of the detector current. Since 


lim Ii{t) 

t^oo 


j{st) _ 

-'1 ~ 


lim s /i(s) = So 

s—>0 


Ad A $ 1 ( 0 ) 

Si S2 


— So 


Ad ‘hi(O) 


(52) 


it is obvious that a stationary expectation of the detector current exists only when p < 0, i.e. 
when the multiplying assembly is in a subcritical state. The quantity $i(0) characterizes 
the average value of the electrical charge produced in the detector during the registration of 
one neutron. 


3.1. A concrete example for the expectation 


Eqs (51) and (52) show that the expression for the expectation of the detector current 


contains the average current pulse $i(t), generated by the detection of a single neutron 
at time f = 0. In general, the shape of ^i(t) depends on a number of physical processes 
taking place in the detector during the rather complicated processes of charge generation 
and transport. However, in the present work, like in its predecessors, we will only account 
for the fluctuations of the detector current due to the randomness of the arrival times of the 
neutrons. Therefore, similarly to Refs [7] and [1], we will choose a constant value instead 
of the random variable dehned by (|^, i.e. we will use the density function w{x) = 5{x — ae) 
in the formula From this it follows that 


$,(«) = 


(53) 


Based on the shapes of experimentally observed current pulses it appears that the empirical 
expression 

$i(f) = (p(ae,^) = (54) 

is an acceptable approximation, hence it will be used for our illustrative calculations. It is 
seen that plays the role of the decay constant of the detector pulse, and Q is the mean 
value of the charge collected in the case of the detection a single neutron. 


Fig. B illustrates the shape of t 
various values of cxe, whereas Fig. [2 


le time dependence of a single average current pulse for 
shows the stationary detector current during a time 


interval for an aggregate of several pulses. Obviously, the expectation and the variance of 
the are constant. 

In order to evaluate (46) with the <h(t) of (53), one needs the Laplace transform of *hi(t), 
which is obtained as 


4>i(s) = 


e-"‘<hi(f) dt = 


(y.p_ 


s + a. 


Q. 


Hence one arrives at 


ii(s |n(0) = 1) = 


Ad {s + A) 

(s + Si) (s + S 2 ) + Oi 


Q, 


(55) 


(56) 
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Figure 1: Time dependence of a single average current pulse 



Figure 2: Stationary detector current 


Vi 


(st) 


in a given time interval 


where si and S 2 are defined by (|47|) and (|48|), respectively. The inverse Laplace transform of 

S2 - A 


(56) is obtained in the form 
ii{t |n(0) = 1) = XaalQ 


Si — A 


e-sit _ 


o—S2t 




.(si - S2) (Sl - fte) (Sl - S2) (S2 - Cte)^ 

(2 — 52^) Oe A — Ctg t — S 2 A + Ctg (— 1 + S 2 t + A t) 


+ 


-ctet 


(57) 


(Si Cte) (S2 Ctg) 

Si {cte t (cKg — A) — A + S 2 [1 + t (—Cte + A)]} 1 

(Si — CTe) (S2 — Cte) J 

Fig. I displays ( [5^ , showing the time dependence of the expectation of the detector 
signal dne to the nentron chains indnced in a snbcritical assembly by one single starting 
nentron, for three different reactivities. 

The time evolntion of the detector cnrrent from time t = 0, at which time the external 
sonrce was switched on in a system that previously did not contain any neutrons, is given 
by 


Ii{t) = So / ii{t' |n(0) = 1) dt'. 


(58) 


Fig. I^shows the expected effect of the delayed neutrons, forming an intermediate plateau- 
like part in the curves approaching the stationary mean value of the detector current. In 
order to calculate the the stationary expectation of the detector current, i.e. the quantity 

= lim 


t^OO 


(59) 
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Figure 3: Time dependence of the expectation of the detector signal due to the reaction chain induced in a 
subcritical assembly by one single starting neutron, for three different reactivities 



Figure 4: Approach to the stationary value of the expectation of the detector current in a subcritical assembly 
with a source switched on at t = 0, for three different reactivities 


one can use the Tauber theorem as in (52). Taking into account that $ 1 ( 0 ) = Q, one has 

(60) 


= So —^ Q, where p = l3 — aA<0. 


-p/A 


4. Variance of the detector current 


By using (30) and (31), one can write the variance of the detector current for t > 0 in 
the form: 


{77} = K2{t) = hit) = So - 


1 

to 

d'^gioj, t' n(0) = 1) 

Jo 

du^ 


dt', (61) 


- tj=0 


where 


d'^gicOyt' |n(0) = 1) 
duh 


= hit' |n(0) = 1) 


(62) 
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is the variance of the detector current generated by a single starting neutron. An equation 
for i 2 {t |n.(0) = 1) can be derived from (28). After a simple algebra one obtains 


kit |n(0) = 1) = 


Xd / ^2it') dt'+ Xf / 



d‘^qi(^,t') 


- CJ =0 


where 


and 


fiV 

d'^hioj, t') 

vv 

doj^ 


dt\ (63) 


(64) 


q(u, (') = [g(u, t' |n(0) = 1)] [g{u, t' |c(0) = 1)]. 


J aj=0 




(65) 


In the hrst step one writes 

^ ^ d‘^qi^,t') 

doj^ 


iV f [gii^.f |n(0) = 1)] 


a ;=0 



dur 


+ 


ui=0 


2dqsd.(nn(0) = i)..(nc(0 ) = d + ())' javyw^^two)^ 


cj =0 


where 


iX^ j [giu,t'\niO) = 1)] 


= q? [hit' |n(0) = 1)]^ + q'f’ hit! |n(0) = 1), 


2 , hp)(+' 


doj"^ 


aj=0 


and 



= qf^ [hit' |c(0) = 1)]^ + qf^ hit' |c(0) = 1). 


cj=0 




= q:? [hit' |n(0) = 1)]^ + q[^^ hit' |n(0) = 1) + 


aj=0 


2q[^^ q[‘^Uiit'\niO) = l)A(t'|c(0) = 1) + 

q? [hit' |c(0) = 1)]^ + hit' |c(0) = 1). 
Before calculating the function hit |?^(0) = 1), the parameters 


( 66 ) 


(p) (p) (d-) (d) 

q{\ qi X qi X q^ 
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have to be determined. By using the generating functions q^^\z) and defined by 

formulae (p!^, one obtains 


dz 

d^q^\z) 

dz^ 

dq^^\'. 

dz 

d:^q^^\z) 
dz^ 


J Z = 1 


J Z = 1 


J 2 = 1 


(p) 

= Qi = 


= q2^ =< - 1) >, 


(d) 

= qi = ^d, 


= q? =< M^d - 1 ) > 


J 2 = 1 


In the next step one substitutes (66) into (63), which leads to 

rt 


i 2 {t |n(0) = 1) = A, y $ 2 (t') dt' + Xf j ^ 2 (^' |n(0) = 1) + 

i 2 {t' |c(0) = 1) + q 2 ^ |n(0) = 1)]^ + 

2 q[^^ q[^Ui{t'\n{0) = l)ii(t'|c(0) = 1) + [zi(t'|c(0) = 1)]^| dt'. (67) 


In the further calculations one needs the relations 

rt 


ij{t |c(0) = 1) = A / e ij{t' |n(0) = 1) dt', j = 1, 2, 


( 68 ) 


'0 


which can be obtained from (29). For simpler notation it is useful to introduce the function 

Rii') = q? Pi(^' 1^(0) = 1 )]^ + 

2g{^^ q[‘^^ ii{t' |n(0) = 1) ii{t' |c(0) = 1) + qi"^^ [ii{t' |c(0) = 1)]^ (69) 

into the integral equation ( [^ . Eq. ( [^ , can be solved by Laplace transform. One obtains 

(s + Xr) * 2(5 |n(0) = 0) = 


where 


and 


(p) 

qr 


Xd ^ 2 ('S) + Xf uq * 2(5 1^.(0) — 1) + Aj Ufi * 2(5 |c(0) — 1) + Aj R{s), 


i 2 {s |c(0) = 1) = ^-r ^ 2(5 |n(0) = 1) 
s + A 


(70) 


(71) 


ms) = / e-“ R{t) dt = 

Jo 

poo 

[ii(t |n(0) = 1 )]^ dt + 2 q^f'^ q['^^ / Zi(t |n(0) = 1) ii(t |c(0) = 1) (it+ 
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(72) 


(d) 


e [ii{t |c(0) = 1)]^ dt. 

Jo 

Accounting for ( [7l| , after elementary calculations one arrives at 

A 


S + \r - XfUo - XfUd 


s -|- A 


where 


i.e. 


S + Xr - Xf Uo - Xf Ud 
i 2 {s |n(0) = 0) = 


^ 2(5 |n(0) = 0) = Ad $ 2 ( 5 ) + XfR{s), 


X (-s + Si) (s + S 2 ) 


s -|- A 
s “1“ A 


(s + Si) (s + S 2 ) - 


s + A ’ 
Xd *^ 2 ( 5 ) + Xf R{s) 


(73) 

(74) 

(75) 


By using (61) and (62), the Laplace transform of the variance of the detector current Rit) 
is obtained as 

7 ( \ _ (5 + A) 

R\s) — So 


s (s + Si) (s + S 2 ) - 


Ad $ 2 ( 5 ) + XfR{s) 


(76) 


The calculation of the inverse Laplace transform of Ris) is rather complicated task. It is 
easier to calculate it by symbolic manipulation codes. We have used Mathematica m by 
Wolfram for solving the present problem. 

4.I. Concrete example for the variance 


The second moment of rjsit) with the selected particular detector pulse shape (54) has 


the form 


4>2(t) = [^Rt)f = ate Q\ 

(77) 

whose Laplace transform is 


(78) 

Since the expression of 

R{s) = 

R^^\s) + 2 qf ''+ q^^'^ R^'^\s) 

(79) 


is prohibitively long, we do not reproduce it here in print. 

In order to illustrate the effect of prompt neutrons just after the switching on of the 
neutron source, the time dependence of R{t) for small values of t was calculated. Fig. 
shows the time dependence of the variance R{t) of the detector current in a subcritical 
assembly just after the neutron source was switched on, for three different reactivities. From 
this hgure it would appear as if the system, as monitored by the detector signal, reached the 
stationary state rather fast. An inspection of the long-time behaviour of the system shows, 
however, that this is note the case, and one has to follow up the behavior of the function 
R{f) during a much longer period to arrive at the stationary variance of the detector current. 

Fig. @ shows the approach to the stationary value of the variance of the detector 
current for three different reactivities during a longer time period. It is seen how the presence 
of delayed neutrons extend signihcantly the time needed to reach stationarity. 
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Figure 5: Time dependence of the variance of the detector current l 2 it) in a subcritical assembly just after 
the neutron source was switched on, at three different reactivities 



Figure 6: Approach to the stationary value of the variance of the detector current in a subcritical assembly 
with a source switched on at t = 0, for three different reactivities 


In practice one needs the dependence of the stationary variance of the detector cnrrent 
on the reactivity p. With the help of the Tauber theorem, one hnds that 


D' 


\d OLp A 


^ jlso ^ lim 5 / 2 ( 5 ) = 5 o , 

^ ' s -!>0 45i 52 


0"C(p), 


(80) 


The explicit form of the function C(p) is rather lengthy, therefore it is given in the Appendix. 


The reason for separating out the multiplying factors in (80) from the function C(p) will be 
clear in the forthcoming discussion, where a comparison with the results of the traditional 
Campbell theorem will be shown. For an illustration, the dependence of I 2 ^ on the reactivity 
p is shown in Fig. 

Fig. shows the sensitivity of to the variation of the detector pulse decay constant 
tte- One hnds that a larger cte brings about a larger stationary variance of the detector 
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Figure 7: Dependence of the stationary variance of the detector current on the reactivity p of a 
subcritical assembly 




Reactivity p 


Figure 8: Dependence of the stationary variance of the detector current I 2 on the reactivity p of a 
subcritical assembly at two decay constants 


current. 

5. Discussion 

In possession of the result for the variance of the detector signal for the case of detection 
in a multiplying medium, it is possible to compare it with the value obtained from the 
application of the traditional Campbell formula. Such a comparison was already made in 
our previous work [T], for the case where only prompt neutrons were assumed in the hssion 
chain. In addition to performing the same analysis by accounting for the delayed neutrons, 
also some further aspects will be discussed, which were not analysed in the previous work. 

As mentioned in [T], for a correct comparison, one has to account for the fact that the 
traditional Campbell formulae are expressed in terms of the intensity of the detection events 
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in the detector (which will be denoted here as Sq), whereas in the present formulae the 
intensity sq of the injection of neutrons from an extraneous source appears. The asterisk 
here is to indicate that the two intensities are intrinsically different. In the continuation, in 
order to distinguish between the traditional formulae (independent events) and the present 
ones (non-independent detection events), the former will be denoted by an asterisk. 

Since in a subcritical medium with reactivity p and an extraneous neutron source with 
intensity sq, the stationary neutron density is given a^ = soA/(—p), the detection 
intensity Sq will be equal to 

So = ^0— (81) 

For a correct comparison, the traditional Campbell formulae need to be used with a detection 
intensity Sq above. 

The stationary variance of a detector signal with independent incoming events, with 


the detector response function given by Eq. (54), was already calculated in Ref. [7j with the 
result 

■r»2 / (®*) 1 ^ * /n2 ^ ^dAo.^ 2 ^dOie^ ^2 

D i r = = 7 S 0 - Q = So - - Q . (82) 

l J A 4 —p 4si S 2 

Here, in order to facilitate the comparison with the formula obtained in this paper for non- 
independent incoming events, Eq. (80), in the last equality we used the identity S 1 S 2 = 
—Xp/A, which can easily be obtained from Eqs (47) and (48). With Sg properly expressed 


in terms of Sg, the bias of the traditional formula, when using it for the evaluation of 
measurements made in a subcritical or a critical core where the primary detection events 
are not independent, can be expressed by the ratio of the correct variance obtained in the 


present results, to the variance of the traditional Campbell method. By using (80) and (82) 
one obtains 

D2 


DMr; 


4*0 


= C(rt 


(83) 


where the function C(p), giving the bias of the traditional formula, is given in the Appendix. 
Figure]^ shows the dependence of the bias factor C(p) on p. The sensitivity of the bias 


factor to the detector signal decay constant cxe can be seen in Figure 10 


The dependence of C(p) on p, although algebraically much more complicated, follows the 
same tendency as in the previous case, i.e. without delayed neutrons in the multiplication 
process (Ref. [1], Fig. 7). It is seen that, as the system approaches criticality ((—p) —t 0), 
the bias factor diverges. This behaviour can also be readily derived from Eq. (| 88 |), by noting 


from (47) that 


lim Si = A ^ 

{-p)^o A 


and 


S2 —t 


A 


AA + /3 


p if — p — 0 . 


(84) 


Using also the identity 


S1S2 = 


' A 


^For a recent note on a general misconception regarding the derivation of this formula, see Ref. 
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Figure 9: Dependence of the factor C(p) for the stationary variance of the detector current on the reactivity 



Figure 10: Sensitivity of the factor C(p) to the decay constant Oe 


it is readily seen that for vanishing reactivities, C(p) diverges as (—p)“h This is becanse 
the stationary variance of the nentron number in a critical reactor diverges, and the correct 
formula for the variance of the detector signal mirrors this fact. 


From the practical point of view, Eqs (80) and (81) show that in measurements in a 


subcritical core, the variance of the detector signal is still proportional to the detection 
intensity, which in its turn is proportional to the neutron flux, although with a proportionality 
factor which is not given properly by the traditional Campbell theorem. In addition, the 
bias of the traditional Campbell formula depends also on the level of subcriticality. However, 
from the quantitative point of view. Fig. shows that even for moderate subcriticalities (i.e. 
close to critical, such as p = 0.002 —0.3$, corresponding to kg/f = 0.998), the bias factor 

is still quite close to unity. Hence, hssion chambers can be used without problems to measure 
the stationary flux in subcritical systems, as long as the detector was properly calibrated in 
measurements, to obtain the correct proportionality factor. 

The apparent problem that the bias factor diverges for critical systems, does not consti- 
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tute a problem in practice either. In reality, one determines an estimate of the exact value of 
the calibration factor, whose definition is based on an ensemble average, from a time average, 
which is taken over a finite time interval, hence it remains hnite. This is very much in line 
with the way how the auto power spectral densities of time-resolved measurement signals 
are determined from the Fourier transforms of the signal (as opposed to that of its auto¬ 
correlation function) by the help of the Wiener-Khinchin theorem, which also involves the 
appearance of a scaling factor which is inhnite, but whose estimate remains a hnite quantity. 

Interesting insight can be gained by taking the opposite case of deeply subcritical systems 
with p —)■ —oo. Fig. |9]mdicates that in that case the bias factor tends to a constant. A 
simple analysis of Eq. ^88| shows that, since 

Si —)■ — Y if P —t —oo; and lim S 2 = A (85) 

A p ^-00 

one has 

lim C(p) = 1 (86) 

p—^—oo 

This is a completely logical result, which expresses the fact that in a deeply subcritical system 
the effect of branching diminishes. Hence, in the limit, the individual detection events will 
become independent, and the traditional formula and the one accounting for the correlations, 
give exactly the same result. This shows a nice and reassuring link between the traditional 
Campbelling theory and the more involved case where the non-independent character of the 
incoming primary events is accounted for. It is worth mentioning that the same agreement 
between the traditional and the extended theory is found in our previous paper [T], i.e. that 
the bias factor converges to unity when p —)■ — cx). This can be obtained analytically from 
Eq. (60) of Ref. P, which readily yields 

lim C(a) = 1 (87) 

However this fact was failed to be mentioned in p. 

6. Conclusions 

The previously introduced formalism for the calculation of the statistics of the signal 
of a hssion chamber, detecting neutrons in a multiplying medium and hence experiencing 
non-independent detection events, was extended by accounting for delayed neutrons. The 
variance of the detector signal was derived and explicitly calculated with the assumption of 
a plausible detector response function. A comparison with the variance of the traditional 
formula, given by the Campbell theorem, made it possible to quantify the dependence of the 
bias on the subcriticality of the multiplying system. As expected, the deviation between the 
two cases vanishes in the case of deep subcriticalities, since in a non-multiplying medium 
only the source neutrons will be detected, which are emitted independently from each other. 
When approaching criticality, the effect of the multiplication and hence that of the non- 
independent character of the detections will increasingly dominate, thus the bias of the 
traditional Campbelling technique will increase. In practice, as long as the fission chamber 
is calibrated from measurements, this does not pose a problem. Although, according to the 
results, the calibration factor depends on the system subcriticality, our results show that 
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the bias factor is quantitatively quite close to unity for the regimes in which the planned 
subcritical accelerator driven systems are planned to be operated, hence hssion chambers 
can be used for flux monitoring. For the case of measurements in a critical system, which is 
the most important mode of operation of a hssion chamber, the divergence of the bias, and 
hence that of the corresponding calibration factor is handled by estimating the variance of 
the detector signal from a measurement of hnite time duration. 


Appendix 


The factor C(p) in the stationary variance of the detector current (80) is given by the 
following formula: 

C(p) = 

-;-ry- ■®2 ('S 2 + c^e)^ + sf S 2 (s 2 + (s 2 + 2ae) + 

S 1 S 2 (S 1 + S 2 ) (si + ae) (s 2 + ae) 

2 Arf A/ (2 + q^^^ + ttg (§2 + aef + 

2sl S2 tte + 5 si si al + 4 sl sjal + 2 Xd Xf (2 qf^ qf'* + q^f^ + qf^^ Si+ 

sj S2 al + XdXf qi^^ al + XdXfX'^ (2 q^^^ + q^^^ + ) + 

Si si aj + 2si S 2 Og + 4 Ad A/ A^ si ^2 g|^^ gj"*^ + g^^^ + g^"*^ j «£ + 


( 88 ) 


2XdXf Si S2 2(2 qP q^f^ + qf'‘) A^ + g^^^ (2A^ + a^) 


(d) 


Xp) 


O^eX- 


Si S 2 


^4 I \ \ „(P) ^2 \ \ \ J I (PJ (aj I (pj I l( 

«£ + -^d <i 2 ttg + Ad A/ A ( g^ qi + q 2 + q 2 


Sp) Jd) 


Xp) I 


where Si and S 2 are dehned by (47) and (48), respectively. It is worth noting that the 
formula depends, among others, on the detector characteristics (the pulse shape and the 
corresponding decay constant cce of the detector), hence it cannot be considered as universal. 
However, the qualitative monotonic behaviour, as well as the asymptotic properties for the 
cases lim—p —>■ 0 and limp —)■ —00 are universal, and do not depend on the detector 
properties. It is only the speed of the convergence which is dependent on the detector 
characteristics. 
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